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Statistical translation invariance protects a topological insulator from interactions 
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We investigate the effect of interactions on the stability of a disordered, two-dimensional topolog¬ 
ical insulator realized as an array of nanowires or chains of magnetic atoms on a superconducting 
substrate. The Majorana zero-energy modes present at the ends of the wires overlap, forming a 
dispersive edge mode with thermal conductance determined by the central charge c of the low- 
energy effective field theory of the edge. We show numerically that, in the presence of disorder, the 
c = 1/2 Majorana edge mode remains delocalized up to extremely strong attractive interactions, 
while repulsive interactions drive a transition to a c = 3/2 edge phase localized by disorder. The 
absence of localization for strong attractive interactions is explained by a self-duality symmetry 
of the statistical ensemble of disorder configurations and of the edge interactions, originating from 
translation invariance on the length scale of the underlying mesoscopic array. 

PACS numbers: 03.65.Vf, 71.20.-b, 73.20.Fz 


Introduction — In the very first topological insulator 
ever discovered, the quantum Hall insulator [T|, interac¬ 
tions have a dramatic effect by changing the quantum of 
conductance from the value e 2 /h to a fraction of it [2]. 
We know now that the quantum Hall insulator is but one 
entry in a periodic table of topological states of matter, 
including both insulators and superconductors (and com¬ 
monly referred to as topological insulators or TIs) [3-6]. 
These materials are all characterized by a gapped bulk 
with gapless surface or edge excitations, protected from 
localization by disorder due the existence of a topological 
invariant. 

Because a topological invariant is a property of the 
single-particle Hamiltonian, it is a challenge to classify 
TIs in the presence of interactions hhh] Besides the 
question of whether interactions will localize the surface 
states, one can also ask whether interactions may lead to 
distinct topological phases — as they do in the fractional 
quantum Hall effect. Here we address these two questions 
for the superconducting counterpart of the quantum Hall 
insulator, a two-dimensional (2D) superconductor with 
chiral p-wave pairing. 

The superconducting analogue of the electronic quan¬ 
tum Hall effect, the thermal quantum Hall effect, refers 
to a heat current carried by Majorana modes propa¬ 
gating along the edge of a topological superconductor 
[Mm. In the absence of backscattering along the 
edge, the thermal conductance equals G t hermai = cGo, 
with Go = 7T 2 k^T/3h the thermal conductance quantum 
for free electrons. The coefficient c, called the central 
charge, governs the stress-energy tensor of the conformal 
field theory associated with the low-energy edge modes 
[MUIMH]. While c = 1 for both the integer and frac¬ 
tional quantum Hall effects, a Majorana edge has c = 1/2 
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FIG. 1. Array of nanowires on a superconducting substrate, 
with a delocalized Majorana edge mode composed out of cou¬ 
pled zero-modes localized at the end points. 


— at least in the absence of interactions m • One of our 
findings is that moderate repulsive interactions between 
the Majorana fermions drive a transition to an extended 
c = 3/2 edge phase. 

Chiral p-wave superconductors may exist naturally 
(Sr 2 Ru 04 is a candidate material [19, 20]), but a ma¬ 
jor recent development is the search for this exotic pair¬ 
ing in semiconductor nanowires [2TH24] and chains of 
magnetic atoms [25] deposited on a superconductor with 
conventional 5-wave pairing. (Alternative proposals in¬ 
clude chains of atoms in optical lattices [26 and chains 
of magnetic vortices on the surface of a 3D TI E3). 
A two-dimensional (2D) array of parallel nanowires, see 
Fig.0 forms an anisotropic topological insulator, called 
a “weak” TI [2811351 because the Majorana mode prop¬ 
agates only along the edges perpendicular to the wires. 
Each nanowire realizes a Kitaev chain [36], with two un¬ 
paired Majorana zero-energy modes at the end points of 
the wire. These zero-modes overlap to form a dispersive 
ID edge mode of Majorana fermions [37ff39] . Backscat¬ 
tering by disorder is not forbidden, yet this “Kitaev edge 
mode” does not localize HflE]. 

Because the effective boundary theory of the Kitaev 
edge is one-dimensional, it is possible to investigate its 
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behavior in great detail with modern numerical tools 
based on matrix product state (MPS) methods (we use 
both variational infinite MPS methods and the density- 
matrix renormalization group — DMRG). The average 
translation symmetry of the mesoscopic array causes a 
statistical translation symmetry of the ensemble of dis¬ 
order configurations and the exact self-dual structure of 
the interacting edge Hamiltonian. This situation is an ex¬ 
ample of a “statistical topological insulator”, protected 
by a symmetry that is present only on average [40 . We 
find that the exact self-dual structure of the interaction 
Hamiltonian protects the gapless Majorana mode up to 
extremely strong attractive interactions. A gapped phase 
does appear for repulsive interactions, but first the Ki- 
taev edge enters a gapless phase with an unusually large 
central charge c = 3/2. 

Interacting Kitaev edge — The end points of each 
nanowire (labeled s = 1,2,...) in the array of Fig. [l] 
form a ID lattice of self-conjugate Majorana operators 
(7s = 7s) ’ governed by the Hamiltonian 

H = -i ^2 a *7s7s+i - k «7s7s+i7s+27s+3- (1) 

S S 

The a s terms describe hopping along the edge and the 
k s terms describe interactions. Only for odd s do the 
four-Majorana terms have the interpretation of density- 
density (Hubbard-type) interactions [27] . The presence 
of the even s = 2,4,... interaction term is dicated by 
the translation symmetry of the mesoscopic array. The 
fact that all four-Majorana terms appear on equal footing 
determines that the edge Hamiltonian is exactly self-dual. 

We perform a Jordan-Wigner transformation to an 
equivalent spin-| representation, by writing the Majo¬ 
rana operators in terms of Pauli matrices: 

k —1 k —1 

72 k = n a ], 72fc-i = Ofc If • ( 2 ) 

3 = 1 3 = 1 

This transformation splits the coupling parameters into 
even, a e s = ct 2 S , ft® = ft 2 and odd, ct° = a^s+i, ft® = 
ft 2 S +i sets. Statistical translation invariance dictates that 
the even and odd sets are indistinguishable in a clean 
system and have the same probability distribution in a 
disordered ensemble. 

The Hamiltonian 0 of the interacting Kitaev edge 
transforms into 

B = - l^ a s a s - a e a s+l 

s s (3) 

+ 2 ^ K S a s a s+l + 2 ^ a s CT s+2 • 
s s 

This spin model is the self-dual anisotropic next-nearest- 
neighbor Ising (ANNNI) model. The standard (/c® = 0) 
ANNNI model [42l 43 has generically a gapped spec¬ 
trum, except for a critical line with c = 0.5 and a float¬ 
ing phase with c = 1. (See Ref. f44l and l45l for a recent 
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FIG. 2. Quantum phase diagram of the interacting Kitaev 
edge as a function of the relative strength £ of interactions 
and hopping along the edge. The central charges associated 
with the gapless phases are indicated. 


study in the context of Majorana zero-modes.) The spe¬ 
cial feature of the Kitaev edge that protects the gapless 
phase is the equivalence of the even and odd coupling 
terms [40l 41] . The ANNNI Hamiltonian § satisfies 
a corresponding self-duality relation [46], which in the 
context of a spin-^ chain would require an artificial fine- 
tuning of parameters. Here the self-duality is inherited 
from the realization of the Kitaev edge Hamiltonian ([!]), 
where it expresses the natural requirement that the trans¬ 
lation of the Majorana operators by one site, i—q s +i, 
should describe the same physical system. Self-duality 
pins the Kitaev edge at a gapless critical point between 
two gapped phases, protecting it from localization by 
disorder or finite but potentially extremely large inter¬ 
actions. 

Phase diagram of the clean edge — In the absence of 
disorder we may set n e s = = k and ct® = ct° = a (> 0 

for definiteness 07 ]). Then the edge is controlled by 
only one dimensionless interaction-strength parameter 
£ = k/ol. 

We examine the phase diagram of the clean Kitaev 
edge using the evoMPS implementation [481150] of the 
MPS time-dependent variational principle and conjugate 
gradient solver in the thermodynamic limit. To deter¬ 
mine the central charge of the gapless (critical) phases 
we divide the infinite chain into half-chains A and B 
and calculate the entanglement entropy S = — Tr In p a 
from the reduced density matrix pA • The scaling 
S = (c/6) ln£ + const, of S with the correlation length £ 
of the slowest decaying correlation function provides an 
estimate of the central charge 0UE2]. A saturation of 
S with £ on the other hand provides strong evidence for 
a gapped phase. In this way we obtain the three phases 
indicated in Fig. [2j Representative numerical data for 
the floating phase are shown in Fig. [3] 

The non-interacting point £ = 0 corresponds to the 
critical phase of the Ising model, with central charge c = 
1/2. We find that this phase persists up to extremely 
large attractive interactions. Our numerics extends up 
to £ = —100, consistent with a very recent independent 
report [53] of a phase transition into a gapped phase for 
£ ~ —250 (see also [60] 5. 

For £ > 0 a second-order phase transition takes place 
at £ ~ 0.3, leading to a new gapless phase characterized 
by incommensurate charge-density waves. We examine 
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FIG. 3. Correlation length dependence of the entanglement 
entropy for different control parameters f in the “floating” 
phase. The straight line fits to S = (c/6) ln£ + const, indicate 
a central charge of c =a 3/2. Bond dimensions used are in 
the range 18 < D < 102. All states were converged up to an 
effective energy gradient norm < 10 -8 . 

the dominant wavevector k corresponding to modulations 
in the connected <j x and a z spin-spin correlation func¬ 
tions (see Fig.[4|. In the Ising phase (£ < 0.3) these cor¬ 
relation functions have diverging correlation length and 
no modulations, so k = 0. At £ ~ 0.3 the system enters a 
“floating” phase in which the wavevector varies continu¬ 
ously. Unlike the floating phase of the standard ANNNI 
model, the floating phase of the Kitaev edge has has an 
unusually large central charge of c = 3/2, see Fig. |3j. 

For sufficiently large £ > 5, the wavevector locks to a 
commensurate value: k = 7 t /2 for a z and k = tt for a x . A 
finite discontinuity in the second derivative of the energy 
accompanies the locking of the wavevector. We interpret 
the locking of the wavevector as a transition into a phase 
characterized by a commensurate charge-density wave, 
akin to the gapped anti-phase of the standard ANNNI 
model. The value of ( where the locking takes place de¬ 
pends significantly on the size of our numerical simulation 
(more precisely, on the bond dimension D of the infinite 
MPS), and therefore the transition point at £ « 5 is a 
D oo extrapolation, as shown in Fig. [4j 

Effect of disorder — We model disorder in the Kitaev 
edge by shifting the nearest-neighbor hopping terms 
and a e s by a random amount in the range [—5,5], drawn 



FIG. 4. Left panel: Locking of the wavevector corresponding 
to <j x modulations as a function of the control parameter £ 
for different bond dimensions D. Right panel: Extrapolation 
of the locking point, to infinite bond dimension. 


independently and uniformly for each lattice site. Be¬ 
cause <a° and a e s are statistically equivalent, the trans¬ 
lation invariance on long length scales is not broken by 
disorder. The disorder-averaged entanglement entropy S 
of a delocalized ID system with open boundaries (divided 
into segments of length x and L — x) is given by [541456] 

S(x) = \cln [(2L/tt) sin (ttx/L)] + constant. (4) 

This formula is an adaptation of the well-known clean 
case (SUET], with c an effective central charge instead 
of the usual central charge c associated with translation- 
invariant systems. 

Previous work [54J [56] has shown that the addition of 
disorder to a critical Ising phase (c = 1/2) drives the spin 
system to a gapless phase of random spin-singlets, each 
contributing In 2 to the entanglement entropy. Charac¬ 
teristic signatures of this phase are a) an effective central 
charge c = yin2 « 0.347, and b) the appearance of a 
peak at In 2 in the probability distribution of S due to 
the singlet contribution. In Fig. [5] we search for these 
signatures, both in the c = 1/2 Ising phase and in the 
c = 3/2 floating phase, using an MPS implementation of 
the DMRG method j58[ |59| . 

For the Kitaev edge at ( = 0 we find that the effec¬ 
tive central charge converges quickly to c = yin2, as 
expected (green data points in Fig. [5^i) . For small values 
of ( < 0, we also observe the same behavior (blue data 
points in Fig. [5^i) very clearly. The fact that the Kitaev 
edge remains delocalized for small £ < 0, or equivalently, 
that the self-dual ANNNI model remains in the random 
spin-singlet phase, is further confirmed by the narrow 
peak developing at In 2 in the half-chain entropy distri¬ 
bution, see Fig. Ih Finite-size corrections become more 
and more significant with increasing — £, making it dif¬ 
ficult to reach a good convergence for £ = —1. The red 
data points in Fig. [5^t give a value of c which we believe 
has not yet fully converged for L = 300, but still seems 
consitent with a delocalized edge. 

Simulations in the floating phase (0.3 < ( < 5.0) are 
computationally more expensive since the entanglement 
entropy, and therefore the required MPS bond dimension, 
is larger in this phase (keeping L fixed) due to the unusu¬ 
ally large central charge c = 1.5. Moreover, in order to 
avoid incommensurate spin-spin correlations, it becomes 
necessary to select values of n which give a family of sys¬ 
tem sizes L commensurate with the ground-state wave 
vector. Figure [5]: shows that the disorder-averaged en¬ 
tanglement entropy saturates as the middle of the chain 
is approached for the largest system sizes available, in¬ 
dicating a gapped phase, albeit with a large value of S. 
The probability distribution of Fig. it also shows a fi¬ 
nite weight of the distribution for vanishingly small val¬ 
ues of the average entropy, which is another indication 
of a gapped phase. We conclude that the floating phase 
is localized by disorder, most likely due to the pinning 
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FIG. 5. Top panels: Effect of disorder on the Ising phase 
at a = 1 with 5 = 0.5. (a) Scaling of the average entan¬ 

glement entropy for attractive interactions n = 0,—0.1,—1 
at L — 100,200,300 respectively, according to Eq. Q. The 
effective central charge is obtained from the slope of linear 
fits (solid lines) in the limit x —» L/2 as c/6. Data has been 
shifted vertically for clarity, (b) Normalized probability distri¬ 
bution of the half-chain entanglement entropy for n = —0.1, 
showing a developing narrow singlet peak at In 2. Bottom 
panels: Effect of disorder on the floating phase at a = 1, 
S m 0.5, and repulsive interactions n = 0.5. (c) The average 
entanglement entropy for L = 32, 64, 96 saturates for large 
values of the scaling function. Lines are guides to the eye and 
the data has been shifted vertically for clarity, (d) Normal¬ 
ized probability distribution of the half-chain entanglement 
entropy for L = 64, showing evidence for states with entan¬ 
glement entropy S 0. Up to ~ 8000 disorder realizations 
were employed for each simulation. 

of incommensurate charge-density waves by the random 
spatial fluctuations of the disorder potential. 

Summary — We have investigated a strongly- 
interacting topological insulator stabilized by a symme¬ 
try that is broken locally but restored on average. Such 
a statistical topological insulator may be realized at the 
edge of an anisotropic p-wave superconductor. Of par¬ 
ticular interest is the possibility of realizing the inter¬ 
acting Majorana edge mode studied here starting from 
a two-dimensional array of vortex lines [27] . since this 
could produce larger interaction strengths than propos¬ 
als involving semiconducting nanowires or atomic chains 


proximity coupled to superconductors. Unlike typical 
one-dimensional models, the effective edge theory of the 
system remains critical, even for large attractive inter¬ 
action strength and/or disorder strength. This behavior 
can be traced back to the average translation symmetry 
of the two-dimensional bulk, which imposes an average 
self-duality on the strongly-interacting Kitaev edge. We 
hope that our work will motivate the search for other 
strongly interacting topological phases in which average 
symmetries of the lattice lead to boundaries which re¬ 
main delocalized in the presence of disorder. 

Acknowledgements — We thank A. R. Akhmerov, E. 
Altman, P. Fendley, M. Maksymenko, Y. Nakata, Z. 
Nussinov, J. Ruhman, N. Shannon, and T. J. Osborne 
for stimulating discussions. This work is part of the 
DITP consortium, a program of the Netherlands Or¬ 
ganisation for Scientific Research (NWO) that is funded 
by the Dutch Ministry of Education, Culture and Sci¬ 
ence (OCW). ICF thanks the European Research Coun¬ 
cil under the European Union’s Seventh Framework Pro¬ 
gramme (FP7/2007-2013) / ERC Project MUNATOP, 
the US-Israel Binational Science Foundation, and the 
Minerva Foundation for support. AM was supported by 
the ERC Grants QFTCMPS and SIQS and by the cluster 
of excellence EXC 201 Quantum Engineering and Space- 
Time Research. L.S. was supported by an Ali Kaufmann 
fellowship and the US-Israel Bi-National Science Foun¬ 
dation. 


[11 K. v. Klitzing, G. Dorda, and M. Pepper, Phys. Rev. 
Lett. 45, 494 (1980) 

[2] H. L. Stormer, D. C. Tsui, and A. C. Gossard, Rev. 
Mod. Phys. 71, S298 (1999) 

[3] A. Y. Kitaev, AIP Conf. Proc. 1134, 22 (2009) 

[4] A. P. Schnyder, S. Ryu, A. Furusaki, and A. W. W. 
Ludwig, AIP Conf. Proc. 1134, 10 (2009) 

[5] M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 

( 2010 ), 

[6] X.-L. Qi and S.-C. Zhang, Rev. Mod. Phys. 83, 1057 

( 2011 ), 

[7] C.-X. Liu, X.-L. Qi, and S.-C. Zhang, Physica E4 44, 
906 (2012) 

[8] X. Chen, Z.-C. Gu, Z.-X. Liu, and X.-G. Wen, Science 
338, 1604 (2012) 

[9] X. Chen, Z.-C. Gu, Z.-X. Liu, and X.-G. Wen, Phys. 
Rev. B 87, 155114 (2013) 

[10] B. van Heck, E. Cobanera, J. Ulrich, and F. Hassler, 
Phys. Rev. B 89, 165416 (2014) 

[11] T. Senthil, arXiv: 1405.4015 (2014) 

[12] T. Senthil and M. P. A. Fisher, Phys. Rev. B 61, 9690 

( 2000 ), 

[13] N. Read and D. Green, Phys. Rev. B 61, 10267 (2000) 

[14] A. Vishwanath, Phys. Rev. Lett. 87, 217004 (2001) 

15 H. W. J. Blote, J. L. Cardy, and M. P. Nightingale, Phys. 

Rev. Lett. 56, 742 (1986) 

[16] I. Affleck, Phys. Rev. Lett. 56, 746 (1986) 











5 


[17] H. Sumiyoshi and S. Fujimoto, J. Phys. Soc. Jpn. 82, 
023602 (2013) 

[18] B. Bradlyn and N. Read, Phys. Rev. B 91, 125303 (2015), 

[19] A. P. Mackenzie and Y. Maeno, Rev. Mod. Phys. 75, 657 
(2003) | 

[20] C. Kallin, Rep. Prog. Phys. 75, 042501 (2012) 

[21] V. Mourik, K. Zuo, S. M. Frolov, S. R. Plissard, E. P. 
A. M. Bakkers, and L. P. Kouwenhoven, Science 336, 
1003 (2012) 

[22] M. T. Deng, C. L. Yu, G. Y. Huang, M. Larsson, 
P. Caroff, and H. Q. Xu, Nano Lett. 12 , 6414 (2012) 

[23] A. Das, Y. Ronen, Y. Most, Y. Oreg, M. Heiblum, and 

H. Shtrikman, Nat. Phys. 8, 887 (2012) 

[24] L. P. Rokhinson, X. Liu, and J. K. Furdyna, Nat. Phys. 
8,795 (2012) 

[25] S. Nadj-Perge, I. K. Drozdov, J. Li, H. Chen, S. Jeon, 
J. Seo, A. H. MacDonald, B. A. Bernevig, and A. Yaz- 
dani, Science 346, 602 (2014) 

[26] A. Biihler, N. Lang, C. V. Kraus, G. Moller, S. D. 
Huber, and H. P. Biichler, Nat. Commun. 5 (2014), 
10.1038/ncomms5504 

[27] C.-K. Chiu, D. I. Pikulin, and M. Franz, Phys. Rev. B 
91, 165402 (2015) 

[28] K. Nomura, S. Ryu, M. Koshino, C. Mudry, and A. Fu- 
rusaki, Phys. Rev. Lett. 100, 246806 (2008) 

[29] Z. Ringel, Y. E. Kraus, and A. Stern, Phys. Rev. B 86, 
045102 (2012) 

[30] R. S. K. Mong, J. H. Bardarson, and J. E. Moore, Phys. 
Rev. Lett. 108, 076804 (2012) 

[31] L. Fu and C. L. Kane, Phys. Rev. Lett. 109, 246605 

( 2012 ) | 

[32] P. Baireuther, J. M. Edge, I. C. Fulga, C. W. J. 
Beenakker, and J. Tworzydlo, Phys. Rev. B 89, 035410 
(2014), 

[33] Y. Tanaka, Z. Ren, T. Sato, K. Nakayama, S. Souma, 
T. Takahashi, K. Segawa, and Y. Ando, Nat. Phys. 8, 
800 (2012) 

[34] P. Dziawa, B. J. Kowalski, K. Dybko, R. Buczko, 
A. Szczerbakow, M. Szot, E. Lusakowska, T. Balasub- 
ramanian, B. M. Wojek, M. H. Berntsen, O. Tjernberg, 
and T. Story, Nat. Mater. 11 , 1023 (2012) 

[35] S.-Y. Xu, C. Liu, N. Alidoust, M. Neupane, D. Qian, 

I. Belopolski, J. Denlinger, Y. Wang, H. Lin, L. Wray, 
G. Landolt, B. Slomski, J. Dil, A. Marcinkova, E. Mo- 
rosan, Q. Gibson, R. Sankar, F. Chou, R. Cava, A. Ban- 
sil, and M. Hasan, Nat. Comm. 3, 1192 (2012) 

[36] A. Y. Kitaev, Physics-Uspekhi 44, 131 (2001), 

[37] D. Wang, Z. Huang, and C. Wu, Phys. Rev. B 89, 174510 
(2014), 

[38] R. Wakatsuki, M. Ezawa, and N. Nagaosa, Phys. Rev. 
B 89, 174514 (2014) 

[39] I. Seroussi, E. Berg, and Y. Oreg, Phys. Rev. B 89, 
104523 (2014) 

[40] I. C. Fulga, B. van Heck, J. M. Edge, and A. R. 
Akhmerov, Phys. Rev. B 89, 155424 (2014) 

[41] M. Diez, I. C. Fulga, D. I. Pikulin, J. Tworzydlo, and 
C. W. J. Beenakker, New J. Phys. 16, 063049 (2014) 

[42] P. Bak, Rep. Prog. Phys. 45, 587 (1982) 

[43] W. Selke, Phys. Rep. 170, 213 (1988) 

[44] F. Hassler and D. Schuricht, New J. Phys. 14, 125018 

( 2012 ), 

[45] R. Thomale, S. Rachel, and P. Schmitteckert, Phys. Rev. 
B 88, 161103 (2013) 

[46] E. Cobanera, G. Ortiz, and Z. Nussinov, Phys. Rev. Lett. 


104, 020402 (2010) 

[47] The choice a > 0 can be made without loss of generality, 
since its sign can be changed by the gauge transformation 
7 S —y (—l) s 7 s in Major ana language, or cr* —)> — crj and 
<jg (—1 ) s (Jg in spin language. 

[48] J. Haegeman, T. J. Osborne, and F. Verstraete, Phys. 
Rev. B 88, 075133 (2013) 

[49] A. Milsted, J. Haegeman, T. J. Osborne, and F. Ver¬ 
straete, Phys. Rev. B 88, 155116 (2013) 

[50] A. Milsted, J. Haegeman, and T. J. Osborne, Phys. Rev. 
D 88, 085030 (2013) 

[51] P. Calabrese and J. Cardy, J. Stat. Mech. 2004, P06002 
(2004), 

[52] L. Tagliacozzo, T. R. de Oliveira, S. Iblisdir, and J. I. 
Latorre, Phys. Rev. B 78, 024410 (2008) 

[53] A. Rahmani, X. Zhu, M. Franz, and I. Affleck, 
arXiv: 1504.05192 (2015) 

[54] G. Refael and J. E. Moore, Phys. Rev. Lett. 93, 260602 
(2004), 

[55] D. Binosi, G. De Chiara, S. Montangero, and A. Recati, 
Phys. Rev. B 76, 140405 (2007) 

[56] G. Refael and J. E. Moore, J. Phys. A: Math. Theor. 42, 
504010 (2009) 

[57] C. Holzhey, F. Larsen, and F. Wilczek, Nuclear Physics 
B 424, 443 (1994) 

[58] S. R. White, Phys. Rev. Lett. 69, 2863 (1992) 

[59] J. A. Kjall, M. P. Zaletel, R. S. K. Mong, J. H. Bardarson, 
and F. Pollmann, Phys. Rev. B 87, 235106 (2013) . 

[60] See the Supplemental Material. 







6 


Supplemental Material 

We provide more information on the clean Kitaev edge 
in the singular limit ( ± 00 . The dimensionless Hamil¬ 

tonian reads 

h = ±J2 %7s+i7s+27 s +3- (5) 

The two signs are connected by a gauge transformation 
(but only if cr = 0), so we focus on the positive sign. 



Also, for numerical purposes, it is convenient to study 
the Hamiltonian 


H = K e ^2 7s7s+l7s+ 2 7s+3 + 7s7sHrl7s+27s+($) 

s=2k s=2k — l 


as a function of K, e . 



FIG. 6. Ground state energy (left panel), as well as its first and second derivatives (right panel, inset) as a function of K e . We 
have used a constant k q = 1 and bond dimension D — 128 throughout. The first derivative is discontinuous across the phase 
transition (dashed line). 


Figure [6] shows the behavior of the ground state energy, 
as well as its first and second derivatives as K e is varied 
through the self-dual point, keeping k, q = 1. The first 
derivative is discontinuous at the phase transition. To 
further confirm that the transition at the self-dual point 
K e = 1.0 is of first order, we have performed scaling of 
the energy gap as a function of bond dimension. The 
results of Fig. [7] show that the gap remains finite in the 
thermodynamic limit, in agreement with entanglement 
entropy results (see main text for discussion). Note that 
the same behavior was reported in the recent work of 


Ref. |53]. 

Since the Kitaev edge is gapped in the limit ( 00 , it 

is natural to ask whether this phase of infinitely strong 
interactions is adiabatically connected to the gapped an¬ 
tiphase that we find for ( > 5. In order to address this 
question, we investigated the Kitaev edge for > 0. 
We conclude that the antiphase does indeed persist all 
the way to £ 00 , and that there is no transition in the 

singular limit. This situation should be contrasted with 
what happens for £ —> — 00 , as discussed in the main 
text. 
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FIG. 7. Energy difference between the ground state and the 
first excited level as a function of inverse bond dimension, 
1/D, at the self-dual point K e = Ko = 1- The system gap 
scales to a value E\ — Eo ~ 0.092 as the bond dimension is 
extrapolated to infinity, indicating that the self-dual point is 
gapped in the thermodynamic limit. 



